Evidence that FeLoBALs may signify the transition between an 
ultraluminous infrared galaxy and a quasar 

D. Farrah^ 
M. Lacy^ 
R. Priddey^ 
C. Borys^ 
J. Afonso^ 

ABSTRACT 

We present mid/far- infrared photometry of nine FeLoBAL QSOs, taken using 
the Spitzer space telescope. All nine objects are extremely bright in the infrared, 
with rest-frame 1-1000/im luminosities comparable to those of Ultraluminous In- 
frared Galaxies. Furthermore, a significant fraction of the infrared emission from 
many, and possibly all of the sample is likely to arise from star formation, with 
star formation rates of order several hundred solar masses per year. We combine 
these results with previous work to propose that FeLoBALs mark galaxies and 
QSOs in which an extremely luminous starburst is approaching its end, and in 
which a rapidly accreting supermassive black hole is in the last stages of casting 
off its dust cocoon. FeLoBAL signatures in high redshift QSOs and galaxies may 
thus be an efficient way of selecting sources at a critical point in their evolution. 
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Introduction 



First seen in 1967 (lLyndslll967l ). Broad Absorption Line (BAL) QSOs are those objects 
that show broad, deep troughs in their UV and opti cal spectra, arising from resonance hne 



absorption in ga s outflowing with veloc ities of > 0.1c (IWeymann et al.lll99ll : lArav et al.ll2001 
Hall et al.l l2002l : iReichard et al.ll2003l ). BAL QSOs come in three subtypes depending on 



which absorption features are seen. High Ionization BAL QSOs (HiBALs) show absorption in 
Lya, NV A1240, SilV A1394 and CIV A1549, and comprise about 85% of the BAL population. 
Low Ionization BAL QSOs (LoBALs) contain all the absorption features seen in HiBALs, 
and also show absorption in Mgll A2799 and other low ionization species, and comprise 
~ 15% of the BAL population. Finally, a rare class of BAL QSO, in addition to showing all 
the absorption lines seen in LoBALs, also show absorption features arising from metastable 



excit e d levels of iron. These are termed F eLoBAL QSOs (IHazard et al.l 119871 : iBecker et al. 
I997I : branch et aPbooi Ihacv et al.l [2002! ) . 



Efforts to explain the origin of BALs in QSOs have been ongoing since their discovery. 
Broadly, there are two possibilities. The flrst is that BAL QSOs are normal QSOs viewed 
along a particular line of sight; in this case the absorption features arise when an accre- 
tion disk wind encounters a high column d ensity, ionized gas outside the broad line region 
( iMurray fc Chiang] Il998l : iProga et al.l I2OOOI ) . The gas is driven outwards via resonance line 
absorption, but the high column density of the gas shields it from higher energy photons 
that would otherwise completely ionize it. In this case, BA L QSOs are those QSOs viewed 
along a line of sight that coincides with the outflowing gas (lElvisI 120001 ). The second possi- 
bility is that BAL QSOs are youthful objects, still surrounded by gas and dust in which the 
absorption takes place; in this case the BALs do not aris e due to a particular line of sight 
dVoit et ahlEggal : iBecker et al.1ll997l : Iwilliams et"aDll999h . 



There has been signiflcant debate over the years on the best method to flnd young 



QSOs (e.g. ISanders et al.lll988l : iKawakatu et al.ll2006l ). so the idea that some BAL QSOs 



may be such objects is par ticularly intriguing. Most attention has focus ed on the LoBALs 
as candidate young QSOs (ILipari et al.lll994j : ICanalizo fc Stocktonll200ll ). as the differences 
between the line properties of LoBAL QSOs and t hose of ordinary QSOs are hard to explain 
solely in terms of different relative orientations (ISprayberry &: Folta Il992l ). Nevertheless, 
the picture of LoBALs being young QSOs is controversial; for example sub-mm observations 
dLewis et allboosl : Iwillott et aUboOsI : IPriddev et aPboOThsh ow no differences between BAL 



QSOs and ordinary QSOs. Furthermore, IVoit et al.l ( 



19931 ) propose a scenario in which 



LoBALs form via ablation of dust by UV photons in outflows arranged either as a thick disk 
or as an isotropic distribution of clouds; in this case LoBALs do flt within AGN orientation 
schemes, and this scenario is consistent with polarimetric observations (ISchmidt fc Hines 
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19991 : lOde et al.lll999l : iHutsemekers fc Lamvll2000h . 



Recently however, evidence has mounted that FeLoBAL QSOs may be the strong est can- 
didates for being youthful QSOs. Based on UV and optical spectra, iHall et al.l (120021 ) suggest 
that FeLoBAL s are young objects still surrounded by a dust cocoon. Similar conclusions 
are reached by iGregg et al.l (120021 ). who also postulate that FeLoBALs may be associated 
with galaxy interactions. Further evidence comes from observations of the only two known 
FeLoBALs at z < 0.15, both of w hich are associated with Ultraluminous Infrared Galaxies 
(ULIRGs, e.g. iFarrah et al.ll2005l ). Finally, the presence of winds with large spatial extents 
(> lOOpc) in some FeLoBAL QSOs (Ide Kool et al.ll2002l ) provides a plausible mechanism for 
an emerging QSO to directly affect the star formation. Although selection effects may play 
a role, this lends weight to the idea that FeLoBALs and ULIRGs are linked in some way. 

It is plausible therefore that the Fe absorption seen in FeLoBALs arises from iron injected 
into the ISM by an ongoing or recent starburst; marking the FeLoBAL phenomenon as a 
transition stage in a ULIRG when the starburst is at or near its end, and the central QSO is 
starting to throw off its dust cocoon. In this letter we examine the validity of this speculation, 
by observin g a sample of nin e FeLoBAL QSOs using the Multib and Imaging Photom eter 
for Spitzer (IRieke et al.ll2004l ) onboard the Spitzer space telescope ( iWerner et al.ll2004j ). We 
assume a spatially fiat cosmology, with Hq = 70 km s~^ Mpc~^, = 1, and = 0.7. 
Unless otherwise stated, the term 'IR luminosity' refers to the total luminosity integrated 
over 1-1000/im in the rest-frame. 



Observations 



Due to their rarity, assembling a homogenous sample of FeLoBAL QSOs is not trivia l. 
At the time of writing this paper a few hundred FeLoBALs are known (ITrump et al.ll2006l ). 
but when the proposal was being written only 40 or so FeLoBAL QSOs were known over 
the whole sky, most of which lay at z > 1. We therefore aimed to select a sample that 
spanned a relatively narrow redshift range, and gave a reasonable sampling of the properties 
of FeLoBAL QSOs, at the expense of homogeneity, with the broader goal of obtaining a 
qualitative idea of the range of their infrared properties. We imposed an upper redshift cut 
of z = 1.8 and a lower redshift cut of z=1.0. We then randomly selected nine objects. Six 
of these ob j ects w ere found via spectroscopic observations of the Sloan Digital Sky Survey 
(IHall et al.l |2002| ). one (ISO J0056-2738) was di scovered serend ipitously from followup of 
distant clusters that had been surveyed with ISO JDuc et al.lbooi ). one (FBQS J1427+2709) 
was disc overed during spect roscopic followup of radio loud quasar candidates from the FIRST 
survey jBecker et al.l Il997h . and one is the 'archetype' FeLoBAL QSO LBQS 0059-2735 
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( iHazard et al.l 119871 ) . The full sample is listed in Table [H 



Seven of the sample were observed in July 2006, one (SDSS J0338+0056) was observed 
in February 2007, and one (LBQS 0059-2735) was observed in March 2004, using MIPS 
at 24/im, 70//m and 160/im (PIDs 30299 & 82). Following observation, the raw data were 
processed automatically with the MIPS data reduction pipeline at the Spitzer Science Center, 
which performs standard tasks such as image coaddition, sky and dark subtraction, and bias 
removal. We inspected the output frames from this pipeline, and determined that they were 
of sufficient quality for reasonably accurate (~ 10%) photometry of our targets, hence no 
further reduction steps were performed. Photometry was carried out using the 'digiphot' 
package within the Image Reduction Analysis Facility (IRAF) software. We used apertures 
of 2.45, 2.0 and 2.5 pixels at 24;Ltm, 70//m and 160//m to measure the fluxes of our objects, 
resulting in aperture corrections of 1.698, 3.900 and 3.215 respectively. 



Results 



The MIPS fluxes are presented in Tabled! along with any available archival IR photome- 
try. Four objects are detected in all three MIPS bands, though two of these four are detected 
only weakly at 160/Ltm. Three objects are detected at 24/um and 70//m. The remaining two 
objects are detected only at 24/xm. 

In the absence of IR spectroscopy, detailed measurements of the properties of IR- 
luminous AGN or starbursts in our sample are not possible. We can however constrain 
both the total IR luminosities, and the contribution from star formation and/or an AGN, 
by fitting the IR photometry simultaneously with the library of model spect ral energy dis- 
tributions (SEDs) for t he emission from a starburst ( Efstathiou et al.l 2000 ) and an AGN 
( iRowan- Robinson! Il995l ). following the methods described in iFarrah et al.l (120031 ). These 
model libraries span a large number of free parameters (e.g. torus opening angle and line 
of sight for the AGN, burst lifetime and UV opacity for the starbursts) but we here use 
the complete model libraries solely to estimate the likely range in total, starburst and AGN 
luminosities that are consistent with the available data. The results are presented in Table 

El 

The four objects with detections in all three MIPS bands are all extremely luminous, 
with IR luminosities exceeding 10^^'^Lq. In all four objects a starburst is required to explain 
the 160/Ltm emission while remaining consistent with the 24^m and YO/im points (for LBQS 
0059-2735 an AGN model is consistent with the MIPS data, but the detection at 850^m 
requires a starburst component). The predicted starburst luminosities exceed IO^^ ^'Lq in all 
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cases, with inferred star formation rates of several hundred Solar masses per year. The best 
fit SEDs for these four objects are given in Figure [H An example of a pure AGN fit to one 
of these four objects is shown in Figure [21 

The remaining five objects are not detected at 160/im, hence the constraints on their 
luminosities and power sources are cruder. For these five objects we do not present best-fit 
SEDs, but merely summarize the results. Three of these five objects are detected at 24/im 
and 70/um, and all have IR luminosities exceeding lO^^'^L©. The non-detection at 160/im 
combined with the detection at 70/im sets an upper limit on any starburst contribution such 
that an AGN must supply at least some of the 24/im and 70/im emission. Indeed, it is 
possible to explain the 24/im and 70/im emission in all three objects purely with an AGN 
model, though a significant starburst contribution is not ruled out. The final two objects are 
only detected at 24/t m. For one of these (ISO J0056-2738) we also have a 15/im flux from 



ISO (IDuc et al.ll2002l ). allowing us to conclude that the bulk of the IR emission likely arises 
from an AGN. For SDSS J2336-0107 however, we have only a 24/im flux, and hence cannot 
set any meaningful constraints on its IR emission. 



4. Discussion 



Starting with the results in this paper, we examine the conjecture that FeLoBAL QSOs 
represent a specific point in an evolutionary sequence between a ULIRG and a 'classical' 
QSO. We initially only consider the total IR luminosities, listed in Table [2l All of our sam- 



ple are extremely IR-luminous, wit 



(iFarrah et al. 



2003 



Blain et al. 



Lonsdale et al. 



2002 



1 lurn inosities comparable to those of both local ULIRGs 



Borys et al. 



20061) and to high redshift sub-mm bright sources (SMGs, 



20041 ). In most cases the upper limits on the IR lumi- 



e.g. 

nosities exceed IO^^Lq, making the m potentially comparable in luminosity to high redshift 
Hyperluminous Infrared Galaxies flRowan-Robinsonl l2000l : iFarrah et al.l |2002| ) . At time of 
writing there exists no 'classical' QSO sample for which direct comparisons to our sample 
are valid (i.e. one at 1.0 < z < 1.8, matched in optical magnitude and observed with MIPS), 
but it is notable that the IR luminosities of our sampl e subs tantially exceed those of nearly 
all the Palomar-Green QSOs presented by iHaas et al.l (120031 ). 



We next consider the power source behind the IR emission. In 5/9 objects, the SED 
fits predict that the IR emission arises at least in part, and possibly entirely, from an AGN, 
perhaps accompanied by a starburst. In 4/9 objects however, those with detections at longer 
wavelengths, the SED fits demand that a substantial fraction of the total IR emission arises 
from star formation, with implied star formation rates of order a few hundred solar masses 
per year. This implies that FeLoBAL QSOs are preferentially associated with obscured. 
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luminous starbursts that are comparable to the starbursts that are thought to power the 
majority of the IR emission from both local ULlRGs and distant SMGs. There is the caveat 
that our sample is heterogenous, but if we restrict attention to only the six QSOs selected 
from the SDSS, then the fraction of sources with starbursts, at 2/6, remains equivalent 
within the errors, though the small sample sizes render these conclusions tentative at best. 

A potentially more serious caveat however is that the results in Table [2] are dependent 
on the assumption that the SED libraries used in the fitting span the range in SED shapes 
of observed IR-luminous starbursts and AGN. This assumption is a reasonable one, given 
that these same SED libraries give good fits to the SEDs of local ULIRGs, and predict 
starburst and AGN luminosities that are consistent with results from other wavelengths 
( iFarrah et al.ll2003l ). There is, however, one scenario we cannot test for. Although we find 
that a luminous starburst must be present in the four objects detected at 160/im, there do 
exist models for the dust distributions around AGN where the dust is so extended (several 
tens of parsecs) that its outer regions are cold enough to emit significantly at wavelengths 
longward of lOO/zm. Such an extended dust distribution could account for some or all of 
the 160;Ltm emission in these four objects. We cannot formally exclude this scenario, but 
note that interferometric observati ons of local AGN h ave found no evidence supporting such 



extended dust distributions 
possibility to be likely. 



le.g. 



Jaffe et al 



We therefore do not consider this 



Finally, we consider our results together with those from previous work. Prior evidence 
that suggests links between FeLoBAL QSOs and ULIRGs includes; (1) t he only two system s 



arrah et al 



Gregg et al 



2005 ). 



2002|), 



at low redshifts known to contain FeLoBALs are both ULIRGs (e.g. 
(2) FeLoBAL QSOs at high redshifts may be involved in interactions ( 
and (3) large-scale winds in FeLoBAL QSOs m ay provide a mechan ism for the emerging 
AGN to affect star formation in the host galaxy (jde Kool et al.ll2002l ). To these we add (4) 
FeLoBAL QSOs class, extremely IR-luminous, with IR luminosities comparable 

to those of ULIRGs at low and high redshifts, and (5) star formation powers a significant 
fraction of the IR emission in many, and possibly all, FeLoBAL QSOs, with implied star 
formation rates comparable to those inferred for local and high redshift ULIRGs. Overall 
therefore, our results combine with those from previous work to construe strong, albeit 
indirect evidence for evolutionary links between FeLoBAL QSOs and ULIRGs. We therefore 
propose that FeLoBALs mark galaxies and QSOs in which an extremely luminous starburst 
is approaching its end, and in which a rapidly accreting supermassive black hole is in the last 
stages of casting off its dust cocoon. FeLoBAL signatures in high redshift QSOs and galaxies 
may thus be an efficient way of selecting sources at a critical point in their evolution. 
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Table 1. FeLoBAL QSO Sample, and infrared photometry 



Galaxy 


RA (J2000) 


Dec 


z 


rrii 


/l5 


/24 


fro 


/l60 






ISO J005645.1-273816 


00 


56 


45.15 


-27 38 15.6 


1.78 


20.95* 


1.33*^ 


1.56 


<7.0 


<50 






LBQS 0059-2735 


01 


02 


17.02 


-27 19 48.8 


1.59 


17.39* 




7.83 


19.9 


29.0<* 


<168.9 


10.3 


SDSS J033810.85+005617.6 


03 


38 


10.85 


-1-00 56 17.6 


1.63 


18.33 




2.02 


8.4 


55.5<= 






SDSS J115436.60+030006.3 


11 


54 


36.60 


-1-03 00 06.4 


1.46 


17.74 




7.59 


17.9 


<50 






SDSS J121441.42-000137.8 


12 


14 


41.43 


-00 01 37.9 


1.05 


18.77 




4.69 


38.3 


78.9 






FBQS J142703.6+270940 


14 


27 


03.64 


-1-27 09 40.3 


1.17 


18.11 




4.76 


32.6 


68.1 






SDSS J155633.78+351757.3 


15 


56 


33.80 


-1-35 17 58.0 


1.50 


18.01 


8.30^= 


13.93 


24.7 


<50 


<86.1 


<3.6 


SDSS J221511. 93-004549.9 


22 


15 


11.94 


-00 45 49.9 


1.48 


16.49 




10.4 


27.2 


<50 






SDSS J233646.20-010732.6B 


23 


36 


45.10 


-01 07 32.4 


1.29 


21.73 




0.48 


<7.0 


<50 







Note. — All fluxes are quoted in mjy. The 'small' field size was used for all three channels, using the default pixel scale at 
70^m. Exposure times per cycle and number of cycles for all sources except LBQS 0059-2735 were 30s, 10s, 10s and 5, 7 and 
4 respectively, giving on-source exposure times of 2260s, 755s and 85s. For LBQS 0059-2735 the exposure times were shorter; 
3s, 10s and 10s with 1, 1, and 3 cycles respectively. Errors are typically 10% on the 24^m fluxes, 20% on the 70^m fluxes, and 
25% on the 160/^m fluxes. Limits are 3cr. 

HLewis et all ||2003|) 

IDuc et al.l ll2002h 

' IClavell l ll998t) 

'^2a detection 

''2.7a" detection 

*R band magnitude 



Table 2. Infrared Luminosities 



Galaxy 


Lt 


ot 






ISO J0056-2738 


12.48- 


■13.00 


12.40-13.00 


<12.90 


LBQS 0059-2735 


12.95- 


■13.20 


12.80-13.00 


12.44-12.82 


SDSS J0338-(-0056 


12.90- 


■13.40 


<12.40 


12.60-13.40 


SDSS J1154-I-0300 


12.89- 


■13.10 


12.70-12.95 


<12.98 


SDSS J1214-0001 


12.70- 


■12.88 


<12.44 


12.68-12.84 


FBQS J1427-I-2709 


12.81- 


■13.01 


<12.70 


12.70-13.00 


SDSS J1556-I-3517 


13.10- 


■13.30 


12.90-13.23 


<12.50 


SDSS J2215-0045 


13.00- 


■13.40 


12.50-13.40 


<12.80 


SDSS J2336-0107 


11.70- 


■12.90 







Note. — Quoted luminosities are the logarithm of the rest-frame 
1-1000/jm luminosity, in units of solar luminosities (3.826x10^® 
Watts), derived from the SED fits. Limits and ranges are 3a". 
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10 100 1000 10 100 1000 

Fig. 1. — Best-fit observed-frame SEDs for the four objects with detections in all three MIPS 
bands. The best fits have a starburst and AGN component in all four cases, though an AGN 
contribution to the IR emission is required only in LBQS 0059-2735. The fits do not provide 
any constraints on the contribution from a starburst or an AGN in the near/mid-IR, or on 
spectral features (e.g. PAH luminosities). The fits do however illustrate that a starburst is 
required to explain the emission at longer wavelengths, while remaining consistent with the 
24//m and 70/xm fiuxes. Error bars are la while upper limits are 3a. 
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10 100 1000 

A(yU,m) 



Fig. 2. — An example of an observed-frame pure AGN fit to one of the four objects in our 
sample with detections at longer wavelengths, in this case the best possible pure AGN fit to 
SDSS J1214-0001. Here the AGN model differs from that in the lower left panel of Figured! 
the torus is viewed nearly edge-on, leading to a very high inferred obscuration and a strong 
absorption feature at rest-frame 9.7/im. Even with such a model however, the fit misses the 
160/um point by a wide margin. 



